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Abstract. We consider the generation of primordial magnetic fields in a class of bouncing 
models when the electromagnetic action is coupled non-minimally to a scalar field that, say, 
drives the background evolution. For scale factors that have the power law form at very 
early times and non-minimal couplings which are simple powers of the scale factor, one can 
easily show that scale invariant spectra for the magnetic field can arise before the bounce for 
certain values of the indices involved. It will be interesting to examine if these power spectra 
retain their shape after the bounce. However, analytical solutions for the Fourier modes of 
the electromagnetic vector potential across the bounce are difficult to obtain. In this work, 
with the help of a new time variable that we introduce, which we refer to as the e-AA-fold, 
we investigate these scenarios numerically. Imposing the initial conditions on the modes in 
the contracting phase, we numerically evolve the modes across the bounce and evaluate the 
spectra of the electric and magnetic fields at a suitable time after the bounce. As one could 
have intuitively expected, though the complete spectra depend on the details of the bounce, 
we find that, under the original conditions, scale invariant spectra of the magnetic fields do 
arise for wavenumbers much smaller than the scale associated with the bounce. We also 
show that magnetic fields which correspond to observed strengths today can be generated 
for specihc values of the parameters. But, we find that, at the bounce, the backreaction due 
to the electromagnetic modes that have been generated can be significantly large calling into 
question the viability of the model. We briefly discuss the implications of our results. 
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1 Introduction 

Over the last decade or so, a considerable amount of attention has been devoted in the 
literature to study models of the universe wherein, as we go back in time, the scale factor, 
rather than vanish, decreases to a small value, and starts to increase again (see, for instance, 
Refs. [1-14]; for reviews, see Refs. [15, 16]). In other words, in these ‘bouncing’ scenarios, 
the universe goes through a contracting phase, decreases in size to a finite, but non-zero 
value, before it begins to expand. Interestingly, such a behavior for the scale factor helps 
in overcoming the horizon problem, without the need for an inflationary epoch^. One finds 
that, in such scenarios, well motivated, vacuum like, initial conditions can be imposed on the 
perturbations at early times during the contracting phase, as in the conventional inflationary 
scenarios [15, 16]. 

Indeed, there exist (genuine) concerns whether such an assumption about the scale factor 
is valid in a domain where general relativity is expected to fail and quantum gravitational 
effects are supposed to take over. We shall ignore such concerns and refer the reader to the 
literature that discuss such issues [15, 16]. Further, it is well known that matter fields may 
have to violate the null energy condition near the bounce in order to give rise to such a scale 
factor. In this work, we shall not attempt to construct models of matter fields that will give 
rise to a bounce, but shall simply assume such a behavior for the scale factor. 

Various observations point to the presence of coherent magnetic fields in large scale 
cosmological structures as well as in the intergalactic medium (in this context, see Refs. [17, 
18]; for some recent reviews, see Refs. [19-22]). While the field strength in the galaxies and 
clusters is about a few micro-Gauss, an interesting lower bound of a few femto-Gauss on 

^ There will indeed be an epoch of accelerated expansion near the bounce as the universe transits from a 
contracting to an expanding phase. However, the period of acceleration need not last for 60-70 e-folds, as it 
is required in the standard inflationary picture to resolve the horizon problem. 
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the coherent magnetic fields in the intergalactic medium has been obtained [23, 24]. These 
observations cannot be explained by astrophysical process alone and, it seems inevitable that, 
at least on the largest scales, the magnetic fields have a cosmological origin. This has led 
to the construction of various mechanisms for the generation of magnetic helds in the early 
universe. A lot of effort in this direction has focused on the production of magnetic fields 
during the inflationary phase [25-39] . It has long been known that the conformal invariance 
of the electromagnetic field has to be broken if magnetic fields are to be produced in the 
early universe. Else, the expansion of the universe leads to an adiabatic dilution of the 
magnetic field strength to insignificant levels. Often, it is found that, in order to generate 
magnetic fields of sufficient amplitude during inflation, the coupling function which breaks 
the conformal invariance of the electromagnetic action has to grow rapidly at late times. 
Due to this reason, the simplest scenarios of inflationary magnetogenesis have been shown to 
suffer from either the backreaction or the strong coupling problem [31, 40, 41]. 

In this work, we study the generation of primordial magnetic fields in bouncing models. 
We consider a class of scenarios wherein the electromagnetic field is coupled non-minimally 
to a scalar field, which can be possibly driving the bounce. Clearly, it will be interesting to 
examine if scale invariant magnetic fields with sufficient strengths can be generated in such 
models [42, 43]. If one considers scale factors that behave as a power law at early times 
and a coupling function which depends on a power of the scale factor, one can easily argue 
that scale invariant power spectra for the magnetic field can arise before the bounce for 
certain values of the indices involved. An interesting question to ask is if these power spectra 
will retain their shape even after the bounce. Since analytic solutions to the modes of the 
electromagnetic vector potential across the bounce seem difficult to obtain, we investigate the 
problem numerically. We begin by introducing a new time variable which allows for efficient 
numerical integration of the equation of motion governing the vector potential. Imposing 
the initial conditions at sufficiently early times before the bounce, we evolve the Fourier 
modes across the bounce and evaluate the power spectra at a suitable time after the bounce. 
We illustrate that scale invariant magnetic fields are indeed generated under the original 
conditions provided the wavenumbers involved are much smaller than the characteristic scale 
associated with the bounce. 

This paper is organized as follows: In the next section, we shall quickly sketch a few es¬ 
sential aspects of non-minimally coupled electromagnetic fields. We shall discuss the equation 
of motion governing the electromagnetic vector potential, the quantization of the potential 
and the power spectra associated with the electric and magnetic fields. In Sec. 3, we shall 
describe the form of the bounce and the type of non-minimal coupling that we shall consider. 
In Sec. 4, we shall outline the numerical method that we adopt to integrate the equation of 
motion and arrive at the power spectra describing the electric and magnetic fields. We shall 
illustrate that strictly scale invariant magnetic fields arise for a class of couplings and bounc¬ 
ing scenarios. In Sec. 5, we shall examine if magnetic fields that correspond to observable 
strengths today can be generated in the bouncing models. In Sec. 6, we shall discuss the 
issue of backreaction in these situations. Finally, we shall close with Sec. 7 with a summary 
of our results and a brief outlook. 

Before we proceed further, a few words on the conventions and notations that we shall 
adopt are in order. We shall work with natural units such that h = c = 1, and set the 
Planck mass to be Mpj = We shall adopt the metric signature of 

Note that Greek indices shall, in general, denote the spacetime coordinates, modulo A which 
we shall use to denote the polarization of the electromagnetic field. Similarly, Latin indices. 
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apart from k which shall be reserved for denoting the wavenumber, shall represent the spatial 
coordinates. Lastly, an overprime shall denote differentiation with respect to the conformal 
time coordinate. 


2 The non-minimal action, equations of motion and power spectra 

We shall consider a case wherein the electromagnetic field is coupled non-minimally to a 
scalar field ip and is described by the action (see, for instance, Refs. [28, 29, 32]) 



( 2 . 1 ) 


where denotes the electromagnetic field tensor which is given in terms of the vector 
potential as follows: 


( 2 . 2 ) 



The scalar field p could be, for instance, the primary matter field that is driving the back¬ 


ground evolution and J is an arbitrary function of the field. As we had mentioned earlier, in 


order to generate primordial magnetic fields, it is necessary to break the conformal invariance 
of the electromagnetic field. Although there exist a variety of possibilities for breaking the 
conformal invariance, the above non-minimal action is one of the simplest gauge-invariant 
choices. Evidently, it is the coupling function J that is responsible for breaking the conformal 
invariance of the action. 

Variation of the above action leads to the following equation of motion describing the 
evolution of the electromagnetic field: 


= 0 . 


(2.3) 


Consider a (3-|-l)-dimensional, spatially flat, Friedmann-Lemaitre-Robertson-Walker (FLRW) 
universe described by line-element 


ds^ = + o?{t) = a^{'q) {—drf' + dx^^ 


(2.4) 


where t is the cosmic time, a{t) is the scale factor and r/ = / dt/a{t) denotes the conformal 
time coordinate. Let us choose to work in the Coulomb gauge wherein Aq = 0 and SjA* = 0. 


In such a gauge, the equation governing the spatial components of the vector potential is 
found to be 



(2.5) 


In the Coulomb gauge, upon quantization, the vector potential Aj can be Fourier de¬ 
composed as follows (see, for instance, Refs. [29, 32, 40]): 



e*'^- + 6^tA*(r,)e-*’- , (2.6) 


where the modes A^ satisfy the differential equation 



(2.7) 
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In the above Fourier decomposition of the vector potential, the summation over A corresponds 
to the two orthonormal transverse polarizations. The quantities represent the polarization 
vectors, which form a part of the following orthonormal set of basis four vectors: 


The three vectors are unit vectors that are orthogonal to each other and to the wavevec- 
tor k. Hence, they satisfy the conditions 6ij = 1 (with no summation over A) and 
dij A:* = 0. Moreover, the operators and are the annihilation and creation operators 

which satisfy the standard commutation relations, viz. 

[bt bk'] = [bk\ bk>^ = 0 and [b^ b^,^] = 6xx' <5® {k - k') . (2.9) 

If we further define a new variable Ak = J Ak, then the equation (2.7) for Ak simplifies to 

^ Alfc = 0. 


( 2 . 10 ) 


A'l + e- 


r 

j 


One can show that the energy densities associated with the electric and magnetic fields 
can be written in terms of the vector potential Ai and its time and spatial derivatives as 
follows: 


A' A' 

87ra2 ^ 

(2.11) 

^ ffd glm ^ 

(2.12) 


where jo} denotes the spatial components of the FLRW metric. The expectation 

values of the corresponding operators, i.e. and p^, can be evaluated in the vacuum state 
annihilated by the operator b^. It can be shown that the spectral energy densities of the 
magnetic and electric fields are given by [29, 32] 


Kik) 

r^k) 


d(0|PB|0) 
din A: 


d(0|PE|0) 
din A: 


2'k‘^ a^{r]) 


\Ak{v)\^ 


2tx‘^ a^{rj) 


\A'kiri)\^. 


(2.13a) 

(2.13b) 


These spectral energy densities are often referred to as the power spectra for the generated 
magnetic and electric fields respectively. A flat or scale invariant magnetic field spectrum 
corresponds to a constant, i.e. A:-independent, V^{k). 


3 Modeling the bounce and the non-minimal coupling term 


We shall model the bounce by assuming that the scale factor a(r/) behaves as follows: 


a(7/) 


Oq 


1 + {Tihof , 


(3.1) 


where oq is the value of the scale factor at the bounce {i.e. when p = 0), po denotes the 
time scale that determines the duration of the bounce and p > 0. Such a simple choice for 
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the scale factor leads to a symmetric non-singular bounce for which the Hubble parameter 
H = a'/a'^ vanishes at the bounce. In the absence of any detailed modeling, it is natural to 
expect that r^o is determined by the fundamental Planck scale, i.e. rjQ^ ~ Mp^. As we shall 
point out in due course, wavenumbers of cosmological interest are considerably smaller than 
the scale ko = Therefore, as will be clear from our discussion, our essential conclusion, 
viz. that bouncing models can lead to scale invariant spectra for the magnetic field over 
cosmological scales, remains unaffected even if one chooses to works with an that is a 
few orders of magnitude below Mpj. It is useful to note that the above scale factor reduces 
to the simple power law form with 0 ( 77 ) oc rfP at very early times {i.e. when —77 ;§> 770 ). 

We should stress that the condition p > 0 ensures that the universe does not go through 
an accelerated contraction during its early stages {i.e. as 77 —)■ — 00 ). This, in turn, will 
ensure that, at these early times, the modes of electromagnetic vector potential are well 
inside the Hubble radius, thereby allowing us to impose the standard, sub-Hubble, Bunch- 
Davies initial condition. Actually, it is not essential that the modes of the vector potential 
should be inside the Hubble radius in order for us to be able to impose Minkowski-like initial 
conditions during the contracting phase. It would suffice if the conformal invariance of the 
electromagnetic field, which we shall assume to be broken around the bounce, is restored at 
a sufficiently early epoch. In fact, even this is a strong demand and it is not necessary. As 
the equation (2.10) suggests, we can impose well motivated Minkowski-like initial conditions 
at a suitably early time wherein k'^ ^ J"/J. But, we shall require the modes to be inside the 
Hubble radius at early times, if we expect the scalar perturbations too to be seeded during 
the contracting phase. 

Note that the coupling function J is actually assumed to be a function of the scalar 
field 4> that drives the background evolution. As we have discussed, in this work, rather 
than attempt to model the bounce, we shall assume a specific form for the scale factor 
[viz. Eq. (3.1)]. Due to this reason, we shall also assume a particular form for the coupling 
function J. It seems natural to expect that the coupling function will either grow or decay 
away from the bounce. In order to describe such behavior, we shall conveniently express the 
coupling function in terms of the scale factor as follows: 

J{v) = Joa^{v). (3.2) 

Clearly, for re > 0 the minimum of J{r]) is located at the bounce, while it grows away from 
the bounce. In contrast, for re < 0, J{r]) is maximum at the bounce and decays away from 
it. In the context of power law inflation, for the above coupling function, it is known that 
specific choices for the indices lead to scale invariant spectra for the magnetic field [29, 32]. 
As we shall show later, similar choices for p and re lead to scale invariant magnetic fields over 
cosmological scales in the bouncing scenario too. 

4 Numerical analysis and results 

In this section, we shall first outline the procedure that we shall adopt to numerically integrate 
the equation of motion governing the electromagnetic vector potential and then present the 
results we obtain for the class of bouncing scenarios described above. 

4.1 The numerical procedure 

Let us now sketch the numerical procedure that we shall follow to solve the equation of 
motion (2.7) governing the modes of the electromagnetic vector potential. 
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4.1.1 E-AA-folds as the time variable 


Our first challenge is to identify a convenient independent variable for integrating the differ¬ 
ential equation. At first look, this may not seem to be an important issue, but a suitable 
choice can obviously make numerical integration stable, quick and easy to handle. Let us 
explain. It is well known that it is rather inefficient to integrate the equations of motion 
governing fields and perturbations in a FLRW universe in terms of either the cosmic or the 
conformal time coordinates. Because of this reason, for instance, in the inflationary scenario, 
while integrating numerically, one always uses the e-fold as the independent variable that 
represents time (in this context, see, for instance. Refs. [44-48]). Recall that, the conven¬ 
tional e-fold N is defined N = log (a/ao) so that a{N) = gq exp N, where, evidently, = 0 is 
a suitable time at which the scale factor takes the value oq. Due to the exponential function 
involved, a relatively small range in e-folds covers a wide range in time and scale factor. 

However, the function is a monotonically increasing function of N. As a result, while 
e-folds are useful in characterizing expanding universes, it does not seem to be a suitable time 
variable to describe bouncing models^. In a bouncing scenario, it seems to be convenient to 
choose a suitable variable to be zero at the bounce, with negative values corresponding to the 
contracting phase and positive values characterizing the expanding regime. Further, since 
we are focusing on symmetric bounces, it seems natural to demand that the scale factor is 
symmetric in terms of the new independent variable to characterize time. An obvious choice 
for the scale factor seems to be a{J\f) = gq exp (AA^/2), with M being the new time variable 
that we shall consider for integrating the differential equation (2.7). (The factor of two has 
been introduced in the exponential for convenience.) For want of a better name, we shall refer 
to the variable J\f as e-AA-fold since the scale factor grows roughly by the amount between 
M and {M + 1)^. Needless to add, because of the rapidly growing exponential involved, a 
wide range in scale factor can be covered by a relatively short span of e-AA-folds. 


4.1.2 Imposing the initial conditions and evaluating the power spectra 

In terms of the e-AA-fold, the differential equation (2.7) satisfied by can be expressed as 

2 


CW2 


+ 


1 dH 2 dJ J_ 


dAk 

<W 


+ 


kMV 

gH ) 


Ak = 0 . 


(4.1) 


Since the scale factor and the non-minimal coupling function J are specified [cf. Eqs. (3.1) 
and (3.2)], evidently, the coefficients of the above differential equation are straightforward 
to determine. So, if the initial conditions are provided, the differential equation can be inte¬ 
grated numerically using the standard methods. We shall discuss about the initial conditions 
below. Here, we wish to make a clarifying remark regarding arriving at the Hubble parameter 
in terms of the e-AA-fold. Note that the form of the scale factor is specified in terms of the 
conformal time coordinate [cf. Eq. (3.1)]. The corresponding Hubble parameter H can be 
easily evaluated in terms of the conformal time 77 . In order to arrive at the expression for the 
Hubble parameter in terms of e-AA-fold, we shall require the relation between r] and Af. This 
can be obtained from the fact that g{A/) = gqs^ and the expression (3.1) for the scale 
factor. However, caution should be exercised in arriving at the relation. Since the Hubble 

^One can possibly work with functions after the bounce and e~^ prior to the bounce to describe the 
scale factor. But, the resulting cusp at A = 0 is highly undesirable. 

^To be precise, it has to be called e-Af + ^-folds. But, since ~ 1, for convenience, we shall simply refer 
to it as e-Af-folds. 
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parameter is negative during the contracting phase and positive in the expanding regime, 
one has to choose the suitable root for r]{J\f) in order to describe H correctly. For the scale 
factor (3.1), we find that the function r/(AA) is given by 


v{Af) = ±Vo [[a{Af)/ao]^/P-l} 


1/2 


= ±??0 


^MV{2p) _ 


ll/2 


(4.2) 


with the negative root corresponding to the period before the bounce and the positive root 
after. 

Let us now turn to the initial conditions. During inflation, while imposing the standard 
initial conditions on either the scalar or the tensor perturbations, one requires that the 
modes are sufficiently deep inside the Hubble radius. As we have already discussed, such a 
demand is not necessary in the case of the electromagnetic modes. Since the FLRW universe 
is conformally flat, Minkowski-like initial conditions can be imposed on the modes Ak in a 
domain wherein the non-minimal coupling reduces to unity. In fact, even such a requirement is 
not essential and it is sufficient if there exists an early time during the contracting phase when 
^ J” /J'^. The standard initial conditions, viz. that Ak = l/\/^ and A!f. = —i \/k/2, 
can be imposed at such an initial time. Note that these initial conditions transform to the 
following conditions on the variable Ak and its derivative with respect to the e-AA-fold: 


= ^ _L 

" jm 

dAk i Ni [k 1 dJ{J\fi) 1 

'dN ^ H{Mi) J(AA,) V 2 " J2(AA,) dJ\f Tfl’ 


(4.3a) 


(4.3b) 


where A/) denotes the e-AA-fold when the initial conditions are imposed. 

Numerically, we impose the initial conditions at the earliest time when the condition 
= 100 (J"/J) is satisfied. In Fig. 1, we have plotted the quantity J"/J for certain values of 
the parameters involved. Note that, while J”/ J exhibits a single maximum when J is growing 
away from the bounce, it has two maxima when J decays. We should stress here that, in 
the latter case, to impose the initial conditions, we need to carefully choose the earliest 
time wherein the condition = 100 {J"/J) is satisfied. Moreover, it should be pointed out 
that the initial conditions are imposed at different e-AA-folds for different k. Further, since 
the maximum value of J"/ J is of the order of l/??o = ^O) such a condition will always be 
satisfied by modes with wavenumbers such that k < ko. Actually, we find that modes with 
wavenumbers upto k ~ 0(10) ko satisfy the above condition. We evaluate the power spectra 
for modes with wavenumbers up till these values. 

The last point that we need to discuss concerns the evaluation of the power spectra. 
After having imposed the initial conditions at an early time prior to the bounce, we evolve 
these modes across the bounce and compute the power spectra V^{k) and V^{k) in the 
expanding phase. We evaluate the spectra associated with all the modes at the specific time 
when the smallest scale of interest satishes the condition k'^ = 10^ {J" / J). Again, in the cases 
wherein J decays, one needs to be careful in ensuring that the time is the latest time when the 
condition is satisfied. We should mention here that this condition essentially corresponds to 
choosing a domain soon after the bounce when the amplitude of the electromagnetic modes 
are largely constant for modes such that A: <C /cq. 

"^Actually, a similar condition during the inflationary phase would suffice too. But, during inflation, such 
a condition corresponds to the condition for the modes to be inside the Hubble radius. 
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Figure 1. The behavior of the quantity J"/ J has been plotted as a function of J\f for p = 1 and 
n = Z12 (on the left) and n = — 1 (on the right). As we shall see, these choices for the indices p 
and n lead to scale invariant spectra for the magnetic field. The quantity that we have plotted is 
independent of Jo and ag, and we have chosen r/g = e^® = 7.2 x 10^°. We should mention that this 
choice for yg is arbitrary and it has been chosen for illustrative purposes. We shall work with a more 
plausible value of yg in due course. Note that the maximum value of J" j J is roughly of the order of 
I/ 77 Q. Evidently, while J" j J has a single maximum when n = 3/2, it has two maxima when n = — 1. 


4.2 Results 


Before we go on to discuss the numerical results, let us try to understand a few points 
analytically. As we had mentioned earlier, at very early times (i.e. for —77 S> tjq), the scale 
factor (3.1) simplifies to the power law form a(r/) oc rfP. During such times, the non-minimal 
coupling function J also has a power law form and it behaves as J^rj) oc 77 “, where we have 
set a = 2np. In such a case, we have J” j J ~ a{a — l)/rf‘ and it is easy to show that, 
for k ko, the solutions to the modes of the electromagnetic vector potential can be 
expressed in terms of the Bessel functions Ju{x), exactly as it occurs for a similar coupling 
function in power law inflation [29, 32, 40, 43]. One finds that the solutions can be expressed 
in terms of the quantity Ak as follows: 


Ak{r]) = \/-kr] [Ciik) Ja-i/ 2 {-k^) + C 2 {k) J_a+i/ 2 {-kr])] , (4.4) 


where the coefficients Ci{k) and 6 * 2 (A;) are to be fixed by the initial conditions. Upon imposing 
the Bunch-Davies initial conditions as krj ^ — 00 , one obtains that 


Ci{k) 

C2{k) 


^ — i TT 0/2 



Ak 





cos (tt a) ’ 

giTT (a+l)/2 


COS (tt a) 


It can also be shown that 

A'kiv) -yAkip) = ky/-kr] [C'i(/c) J«+i/ 2 (-fe 77 ) - C' 2 (A:) J_„_i/ 2 (-A: 77 )] . 


(4.5) 

(4.6) 


(4.7) 


The spectra of magnetic and electric fields as one approaches the bounce, i.e. as fc 77 —>• 0 , 
can be arrived at from the above expressions for Ak and A'f, — {J'/J)Ak and the asymptotic 


- 8 - 










forms of the Bessel functions. The magnetic field spectrum can be written as [29, 32, 40, 43] 


K{k) = 


s (av-r- 


(4.8) 


where H ~ {2p/aQ rj) {rjo/rj)'^^, while m = a for a < 1/2 and m = 1 — a for a > 1/2, Also, 
the quantity F{m) is given by 

vr 


F{m) = 


2^™+! r2(m + 1/2) cos2(7rm) 


(4.9) 


It is evident that m = —2 leads to a scale invariant spectrum for the magnetic held which 
corresponds to either a = 3orQ; = —2. The corresponding spectrum for the electric held 
can be obtained to be 

where m = l + aifa< —1/2 and m = —a for a > —1/2, while Q{m) is given by 


Q{m) = 


TT 


(4.11) 


2^™-+3 r^(m + 3/2) cos2(7rm) 

Note that when, a = 3 and a = —2, the power spectrum of the electric held behaves as k~‘^ 
and k^, respectively. This implies that, in the bouncing scenario, one can expect these cases 
to lead to scale invariant spectra (over modes such that A; <C Ajq) for the magnetic held before 
the bounce. The question of interest is whether these spectra will retain their shape after 
the bounce. It is for this reason that we shall consider the combination of the indices p and 
n that correspond to either a = 3 or —2. 

It is also possible to analytically understand the behavior of the mode Aj. and its 
time derivative in the case wherein n is positive. Note that, when n > 0, J"/J has a 
maximum at the bounce. In such a case, for k ko, <C J"/ J around the bounce. From 
equation (2.7), it is clear that in a domain where k'^ can be neglected, we have 


A'k + 2 -r ~ 0 , 


J 


which can be integrated to yield 


A'kir]) ^ 




(4.12) 


(4.13) 


J2(?7) ’ 

where 77 * is a time when kf <C J"/ J before the bounce. The above equation can be integrated 
to arrive at 

Akiv) ^ Ak{rjA + A'kiv*) J d?? = AkivA + A'kir] A a^”(T?A J 


V* 


V* 


where we have set the constant of integration to be Ak{r]A- When a = 3, for the scale 
factor (3.1), we can evaluate the above integral to obtain that 


Ak{r]) ~ Ak{riA + A(.(r7* 


a (??*) % 


^2 n 


5 + 3 ( 77 / 770 )^ 

m [1 + ( 77 / 770 ) 2 ]^ 






- 3tan-i(^ 
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X 


. (4.15) 












We should point out that, while the first term is evidently a constant, the second term grows 
rather rapidly during the contracting phase and less so during the expanding regime. We 
should mention that the above solutions for and A'^ are valid around the bounce until 
the condition <C J"/ J is satisfied. In what follows, apart from presenting the numerical 
results for a few different cases of interest, we shall also compare the above analytical results 
with the numerical ones for a = 3. 

We have numerically integrated the differential equation (4.1) using a Fortran code that 
is based on the fifth order Runge-Kutta method [49]. We have also independently checked 
the results we have obtained using Mathematica. In Fig. 2, we have plotted the evolution 
of Ak and its time derivative A'j^ for two widely different modes and certain values for the 
parameters involved. For reasons outlined above, we have worked with indices n and p 
corresponding to a = 3 and —2. It is useful to note from the figures that, when J has a 
minimum at the bounce {i.e. when n is positive), A^ exhibits a maximum, whereas A^ has 
a minimum at the bounce when J has a maximum [i.e. when n is negative). Moreover, it 
is clear from the figures that, for k ko, there exists a wide domain in time over which Aj. 
proves to be a constant. It is over this domain that we shall choose to evaluate the power 
spectra of the magnetic and the electric fields. In the hgure corresponding to p = 1 and 
n = 3/2 (leading to a = 3) and k <C ko, we have also illustrated the analytical results (4.15) 
and (4.13). It is clear that the analytical result for A'j^ matches the numerical result very 
well. The behavior of A^ is also along expected lines. Around the bounce, the analytical 
solution (4.15) we have obtained mimics the exact numerical result extremely well. 

Let us now turn to the evaluation of the power spectra of the magnetic and and electric 
fields generated in the cases of interest. Let us first list out all the parameters that we have. 
These parameters appear in the functions describing the scale factor a and the non-minimal 
coupling J. The parameters that characterize the scale factor are oq, rjo and p. Apart from 
these three, two additional parameters, viz. Jq and n, are required to describe the coupling 
function J. As we have already discussed, k^ = determines the scale associated with the 
bounce. We shall plot the power spectra in terms of k/k^. Clearly, our first goal would be 
to examine if we obtain scale invariant spectra for the magnetic field for any values of the 
parameters. As we have already discussed, for A; <C fco, we expect scale invariant spectra for 
the magnetic field before the bounce when a = 3 or —2. Also, based on the same arguments, 
one can show that, before the bounce, we can expect the power spectrum of the electric field 
behave as k'^~‘^°‘ when a = 3 and as A;®'*'^" when a = —2. With the numerical tools at hand, 
it is interesting to examine if these power spectra retain their shape after the bounce as well. 
In Fig. 3, we have plotted the power spectra V^[k) and V^ik) for a set of cases that lead 
to scale invariant spectra for the magnetic field. It is clear from the figure that the spectra 
indeed retain the shape for small wavenumbers (expected from the analytical arguments) 
as they emerge through the bounce. An interesting point to note from these figures is the 
behavior of the spectra for k > kQ. These modes are hardly affected by the background 
and they retain their original form determined by the initial conditions. As a result, the 
Minkowski-like initial conditions that we have imposed imply that both the magnetic and 
electric fields should behave as k^. This is exactly the behavior that we obtain from the 
numerical results. The fact that the analytically expected results are reproduced indicates 
the robustness of the numerical procedures that have been adopted. 

Note that the dehnitions of the power spectra (2.13) contain an overall factor of J^. 
Further, the differential equation (4.1) satisfied by A^ only involves the ratio (dJ/dW) (1/T). 
Moreover, note that the initial conditions (4.3) contain a factor of J in the denominator and 
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Figure 2. The behavior of the absolute values of Ak (on the left) and its derivative A'/, (on the right) 
has been plotted for the modes fco = Vo^ ~ = 1.389 x 10“^^ (the first and the third rows) and 

k = 10“^° ko (the second and the fourth rows) for the cases wherein n = 3/2 (the top two rows) and 
n = —1 (the bottom two rows). We have set p = 1, oq = 10“^° and Jq = 10"^ in all the cases. It is 
interesting to note that the amplitude of Ak has its maximum value at the bounce when J itself has 
a minimum, while Ak exhibits a minimum when J has a maximum at the bounce. In the absence 
of any detailed modeling of the bounce, we had assumed that = Mp^. So, it is indeed rather 
extreme to consider the mode fco = as scales of cosmological interest correspond to at least 50 
orders of magnitude smaller than kgl However, the evolution proves to be rather simple for much 
smaller wavenumbers. We find that, for fc <C fco, the absolute values ot Ak essentially grows (when n 
is positive) or decays (when n is negative) from its initial value to turn to a constant as one crosses 
the bounce, before it decays or grows with superposed oscillations at very late times. The dashed red 
curves represent the analytical estimates (4.15) and (4.13) and they match the numerical results as 
described in the text. 
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Figure 3. The power spectra of the magnetic (in blue) and the electric (in red) fields for the cases 
wherein (p,n) = (1,3/2) (top left), (p, n) = (2,3/4) (top right), {p,n) = (1,-1) (bottom left) and 
{p,n) = (2,—1/2) (bottom right). Evidently, the top and the bottom rows correspond to a = 3 
and a = —2, respectively. We should stress that we have evaluated these spectra after the bounce. 
We have worked with the same values of rjQ, oq and Jq, viz. po = e^® = 7.2 x 10^°, oq = 10“^° 
and Jq = 10^, in all the four cases. As expected from the analytical arguments, the spectra for 
a given a = 2np [i.e. within each row) have the same shape. However, note that, for a given 
a, the amplitude differs for different p and n. Clearly, the spectrum of the magnetic field is scale 
invariant in all these cases. The power spectra of the electric field are along expected lines, behaving 
as = k~^ when a = 3 and when a = —2 (indicated by the dotted green lines). 

For wavenumbers much smaller than the bounce scale fco, the power in the electric field is smaller 
than the power in the magnetic field when a = —2, i.e. when n is negative. Note that, modes with 
wavenumbers comparable to or greater than the bounce scale fcg a-re hardly affected by the bounce. 
Hence, over these scales, the power spectra of both the electric and magnetic helds are expected to 
behave in a similar manner (in fact, as k^, indicated by the dotted orange lines), which is exactly the 
dependence that we obtain from the numerical results. 
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will hence will involve a factor of 1/Jo- Therefore, the power spectra are independent of Jq. 
(In fact, J also contains the parameter oq. But, it can be absorbed in the overall constant 
Jo leading to the same conclusions as above.) This can also be easily confirmed with the 
numerics. The amplitude of the power spectra are determined by the parameters ao and r/o 
(apart from the indices p and n which also influence the amplitude). Note that, since it is 
the combination k/ao that appears in the differential equation governing Ak, oo simply sets 
the scale. Therefore, the dominant dependence of the power spectra on oo arises due to the 
factor of that appear in the denominators of the power spectra [cf. Eqs. (2.13)]. Hence, 
we can expect the power spectra to behave as For a given set of the other parameters, 
we find numerically that the amplitudes of the spectra indeed behave in such a fashion. 


5 Generating magnetic fields of observable strengths 


We have established that scale invariant spectra of magnetic fields can be generated in bounc¬ 
ing models. Let us now examine if the strengths of these primordial magnetic fields can 
correspond to observable levels today. 

As the issue of reheating in bouncing models remains poorly understood, we shall simply 
assume that, at some stage after the bounce, the universe transits to the radiation dominated 
epoch. We shall also assume that the coupling function J simultaneously reduces to unity. 
Let the bouncing phase end and the radiation dominated epoch begins at the temperature, 
say. Tend) corresponding to the scale factor, say, Oend- Under these conditions, the power 
spectrum of the magnetic field at the epoch corresponding to Tend, say, is related 

to the power spectrum observed today, say, {k), through the relation 


pcnd(^) ^ptoday^^^ ^ ptoday^^) ^ 


where Ttoday is the temperature today. If we assume that Tend — 10^^ GeV (a choice that 
will be explained below) and, since T^oday — 10“^ eV, in order to generate magnetic fields 
of observed strengths today, i.e. Htoday — 10“^® Gauss (in this context, see, for example. 
Refs. [23, 24, 50-53]), the power spectrum V'^'^{k) should be of the order of 


V^^'^ik) = ( 10 - 1 ®)^ 


f 10^2 eV 
VlO-4 eV 


4 

Gauss^ ~ 10^^ Gauss^. 


(5.2) 


We had discussed earlier that the power spectra behave as Uq^. We find that it is indeed 
possible to produce magnetic fields of such large strengths by working with a suitably small 
value of oq. In Fig. 4, we have plotted the power spectra with the required strength for 
certain values of the parameters. A couple of related points require clarification. We had 
mentioned earlier that the spectra are evaluated at a given time when the smallest scale 
of interest satisfies the condition k'^ = 10^ {J"/J). In plotting Fig. 4, we have chosen the 
smallest scale to be ko. For rjQ = 1, p = I and n = —I, we find that this corresponds to 
evaluating the spectra at the e-AA-fold of roughly M = 5.5. This suggests that we can choose 
flend — exp (5.5^/2) Oo — 10 ® no- If the bounce corresponds to the Planck scale Mpj, then 
Tend can be chosen to be Tend = -^pi (oo/oend) = Mpj/IO® ~ 10^^ GeV. It is this choice that 
we have made above. We should add that a more complete calculation relating the strengths 
of the magnetic fields soon after the bounce and the observed strengths today will require a 
good understanding of the transition from the bounce to the epoch of radiation domination. 
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Figure 4. The power spectra with p = 1 and n = —1, corresponding to a = —2 has been plotted for 
a wide range of wavenumbers (with same choice of colors to represent the spectra as in the previous 
figure). We have set rjo = 1, oq = 4 x 10“^® and Jq = in arriving at these results. As we have 
discussed, the amplitude of the spectra are independent of the choice of Jq, which we also observe 
numerically. Clearly, a suitably small value of ag leads to magnetic fields in the early universe that 
correspond to observable strengths today. Note that we have worked with parameters such that 
a = —2 as it results in electric fields of strengths considerably smaller than that of the magnetic fields 
over observationally relevant scales. 



6 The issue of backreaction 

In this section, we shall discuss the issue of backreaction in the bouncing models of our 
interest. In the inflationary context, as we have discussed, scale invariant magnetic fields 
are generated when the coupling function J either grows with the scale factor or decays in 
certain manner. While the former case suffers from the strong coupling problem [31, 40, 41], 
the issue of backreaction becomes important in the latter. In particular, in the latter case, 
the energy density associated with the generated electromagnetic fields rapidly grow with 
time and can dominate the energy density that has been assumed to drive the background 
evolution (see, for instance. Refs. [33, 35]). Such a behavior is untenable and, for the scenario 
to remain viable, the energy densities in the electromagnetic fields that have been generated 
should always remain sub-dominant to the energy density associated with the background. 
Let us now examine if this condition is satisfied in the scenarios that we have considered 
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here. 


Recall that, from the Friedmann equation, we have 


= = 2>Ml 



( 6 . 1 ) 


where /9bg is the energy density that is driving the background evolution. Given the scale 
factor (3.1), the corresponding energy density can be expressed as 


Phg 


12 p'^M^ 
0-0 Vo 


aoj 


-2 (p+l)/p 


( 6 . 2 ) 


The energy density in a specific mode k of the electromagnetic field is given by 

pl^=V,{k)+V^{k). (6.3) 

Evidently, if the backreaction is to be negligible, we require that pbg > P^^ for all modes 
of cosmological interest. Also, this condition should hold true at all times. However, note 
that, since H vanishes at the bounce, pbg does so too. A priori, it should be clear that any 
non-trivial amount of energy density in the electromagnetic fields that have been generated 
will lead to a violation of the required condition. 

Let us nevertheless estimate the energy density associated with the electromagnetic 
modes that have been created. The dependence of on time for a particular mode can be 
easily arrived at from the numerical solutions we have obtained. In Fig. 5, we have plotted 
Pgg for a relatively large scale mode, along with the background energy density pbg- We 
have plotted the results for values of the parameters that we had considered in arriving at 
magnetic fields of observable strengths in Fig. 4. It is clear from Fig. 5 that the energy 
density in the electromagnetic field (of the given mode) is smaller than the energy density 
of the background at early stages of the bounce. However, as one approaches the bounce, 
the energy density of the electromagnetic field grows rather quickly beyond the background 
energy density, leading to a severe violation of the expected condition. Needless to add, this 
issue of backreaction has to be circumvented if bouncing models are considered to be a viable 
scenario for the generation of observable levels of magnetic fields. 


7 Discussion 

In this work, we have studied the generation of primordial magnetic fields in a class of bounc¬ 
ing universes when the electromagnetic field is coupled non-minimally to a scalar field that 
drives the background expansion. We had restricted ourselves to the consideration of sym¬ 
metric non-singular bouncing models that allow initial conditions on the perturbations to 
be imposed at sub-Hubble scales at very early times during the contracting phase of the 
universe. We found that there exists a class of indices describing the non-minimal coupling 
and the scale factor that lead to a nearly scale invariant spectrum for the magnetic field, 
while the corresponding electric field spectrum is sharply scale dependent. We showed that 
certain values of the parameters involved lead to primordial magnetic fields which correspond 
to observable strengths today. However, unfortunately, the backreaction due to the electro¬ 
magnetic fields that have been generated prove to be substantial calling into question the 
viability of the model. 
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Figure 5. The behavior of the energy density in the electric and magnetic fields for the mode 
k = has been plotted (in blue) along with the energy density of the background (in red). 

We have worked with the same values of the various parameters as in the last figure. In the previous 
figure, we had arrived at magnetic fields that correspond to observable strengths today by choosing to 
work with an extremely small value of uq. Such a small value, unfortunately, boosts the energy density 
in the generated electromagnetic fields. The energy density in the electromagnetic fields proves to be 
unreasonably large at the bounce clearly calling into question the viability of the model. 


Although we have not discussed how to obtain the desired bouncing solution, it turns 
out that bouncing scenarios typically require a violation of the null energy condition at the 
bounce. While this is true for spatially flat or open FLRW universes, the necessity to violate 
the null energy condition can be circumvented with a positive spatial curvature. Since the 
backreaction due to the generated electromagnetic fields is an issue at the bounce in spatially 
flat models, an interesting and possible way to avoid it would be to consider models with a 
positive but small spatial curvature. Such a spatial curvature would require a non-vanishing 
energy density at the bounce and hence may aid in overcoming the backreaction problem. 
We are currently exploring such issues. 
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